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Abstract 


Robust control system analysis and design is based on an uncertainty description, called a linear 
fractional transformation (LFT), which separates the uncertain (or varying) part of the system from 
the nominal system. These models are also useful in the design of gain-scheduled control systems 
based on Linear Parameter Varying (LPV) methods. Low-order LFT models are difficult to form 
for problems involving nonlinear parameter variations. This paper presents a numerical 
computational method for constructing an LFT model from a given LPV model. The method is 
developed for multivariate polynomial problems, and uses simple matrix computations to obtain an 
exact low-order LFT representation of the given LPV system without the use of model reduction. 
Although the method is developed for multivariate polynomial problems, multivariate rational 
problems can also be solved using this method by reformulating the rational problem into a 
polynomial form. 


1.0 Introduction 

Formulation of linear fractional transformation (LFT) models of systems involving 
nonlinear parameter variations is of interest for robust control system analysis and design, as well 
as for control of linear parameter varying (LPV) systems. Moreover, the LFT models should be of 
low order for efficient computation during analysis and design. A matrix singular value 
decomposition (svd) approach was presented in 1985 in references [1] and [2] for computing 
LFT’s for problems involving linear parameter variations. However, construction of low-order 
LFT models for problems involving nonlinear parameter dependencies is very difficult, because it 
is equivalent to a multidimensional minimal state-space realization problem for which there is no 
general theory. The approach that has been taken to date for solving nonlinear parameter- 
dependent problems is to successively decompose the system until all components are linear, and 
then to compute an LFT for each linear component based on the result presented in [1] and [2]. 

The LFT’s associated with each system component are then combined using LFT properties to 
form the LFT model of the full system. Model reduction is usually required using this approach, 
because unnecessary repetitions of the varying parameters usually result. A decomposition method 
for LFT modeling of nonlinear parameter-dependent systems was first presented in reference [3], 
and later refined in reference |4]. This latter paper presented a special decomposition approach 
which reduces the number of unnecessary repetitions of the varying parameters, although model 
reduction is still employed to reduce the dimension of the resulting LFT model of the full system. 

The approach presented in this paper is an extension of the computational approach of 
references [1] and [2] for nonlinear parameter-dependent systems, and is based on reference [5]. 
Specifically, the computational approach is developed for multivariate matrix polynomial problems, 
although multivariate rational problems can be solved using this approach by reformulating the 
rational problem to be in a multivariate polynomial form. Reference [6] presents a method for 
doing this. The LFT modeling approach presented in this paper requires no matrix decompositions 
for multivariate polynomial problems, and achieves a low-order LFT model directly - i.e., without 
the use of model reduction. Moreover, the computations are based on simple matrix operations, 
including the svd and solving linear matrix equations. 

2.0 LFT Modeling Problem Definition 

The LFT modeling problem to be addressed in this paper is defined below. It is assumned 
that the problem to be solved is in a multivariate matrix polynomial form. However, as shown in 
reference [6], multivariate rational problems can be reformulated as multivariate polynomial 
problems and solved using this approach. The problem is stated as follows. 
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Given : A linear parameter varying (LPV) model of a nonlinear parameter-dependent system, as 
represented by the following equation 


X' 

= 

A(8) 

B(8)' 

X' 

= S(8) 

X' 

y_ 


C(8) 

D(8) 

u 

u 


6 = [6 1 ,8 2 ,-,6 m ]ER m (2.1b) 


where S(8) has been separated into nominal and varying components, and the varying (or 
uncertain) component, S A (8), has been formulated as an LFT problem given by the following 
equation 


S i (8)-P 21 A(I-P n A)- 1 P 12 .P 21 (I-AP u r 1 AP 12 (2.2) 

in which each element of S A (8) is a multivariate polynomial function of the varying parameters, 8 

Find : A low-order state-space uncertainty model that satisfies equation (2.2) and is characterized 

by the constant matrices P 2b P 12> and P[ | and the uncertainty matrix A(8), as depicted below in 
Figure 1. 
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Figure 1. LFT Model of the Uncertain System 


The P 22 matrix represents the nominal part of the system, and is characterized by the nominal A, 

B, C, and D system matrices. The S A (8) matrix of equation (2.2) is a known matrix of 
multivariate polynomials based on the LPV model for the system. Formulation of this matrix was 
discussed in reference [6], The LFT model equations associated with Figure 1 are given below. 
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(2.3c) 
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where: A(8) = diag [ 6jl nl , 8 2 I n2 , • • • , S m In m ] E R n A x n A ( 2 .4a) 


m 

n A = 2 n i , n i = dim( Ij ) (2.4b) 

i=l 

The LFT modeling problem consists of solving equation (2.2) for P 21 , P ]2 , and P u over some 
low-order A matrix (as defined by equation (2.4)). This is equivalent to a multidimensional 

minimal state-space realization problem over the m varying parameters in 8. Unfortunately, there 
is no existing minimal realization theory for general multidimensional systems (i.e., for m > 3) that 
can be used in solving this problem. In fact, there are no general minimality tests for 
multidimensional systems given a realization. This paper presents a numerical computational 
approach for solving equation (2.2) for P 21 , P 12 , and P u such that the resulting A matrix is of low 
order. These results are summarized in Section 3. 

3.0 Main Results: LFT Model Computation 

As discussed in Section 2, the LFT problem to be solved is given by the following 
equation: 


S a (6) = P 21 (I-A(8)P 11 )- 1 A(8)P 12 , S A (8) E <P "^"cois (3.1) 

The term Sa(8) is a known matrix function of the normalized uncertain parameters in 8, and P 2 i, 

Pl 2 , and Pi i, are the unknown matrix variables to be determined. The dimension of A(8) must 
also be determined in constructing the LFT model such that the resulting dimension is low-order. 

It is assumed that the functional form of the elements of Sa(8) is multivariate polynomial. 

However, as discussed in Section 2, rational problems can also be solved by reformulation of the 
rational problem (see Reference [6]). 

3.1 Numerical LFT Solution Approach 

As can be seen in equation (3.1), solving for the matrices P 2 i, Pi 2 , Pu and A(8) involves 

the inversion of the matrix [I - A(8)Pi j]. For multivariate polynomial problems, this matrix 
inversion can be exactly replaced by a finite series and an associated nilpotency condition. This is 
expressed in the following equations. 

(I-A(8)Pn)-i = I + (A(S)P„) + (A(8)P n )2 + ... +(A(8)Pn) r (3.2) 

[A(8)Pn] r+1 = 0 (3.3) 

Substituting equation (3.2) into equation (3.1) results in the following equation for Sa(8). 

Sa(8) = P 2 iA( 8)P 12 + P 2 i[A(8)P n + (A(8)P n ) 2 + ... (A(8)P u ) r ]A(8)P 12 (3.4) 
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The first term oil the right side of equation (3.4), i.e. P 21 AP 12 , represents the linear uncertain 
components of Sa(8), and the second term adds in the nonlinear terms. For the case of 
multivariate polynomial uncertainties, the nonlinear terms of Sa(8) consist of crossterms of the 8 
parameters and n th -order terms. Thus, the order (r) of the highest term in the series of equation 

(3.4) is determined by the degree of the highest term appearing in Sa( 8), where crossterm degree 
can be defined as follows. 

degree (81 1 82 2 83 3 . . . &i ') = ( It + | 2 + ... + li ) - 1 ; (3.5) 

Then, the exponent r in equation (3.4) can be defined by the following inequality. 

r =s (rij+ri 2 + ... +ri m )-l (3.6) 

where rp is the maximum degree of 8 ; in S A ( 8 ). 

Since the uncertain system matrix, Sa( 8), has as its elements multivariate polynomial 

functions of 8 , it can be easily expanded in a similar manner as the right side of equation (3.4), 
i.e.: 


S a (8 ) = S (8) + S (8) + • • • + S (8) 

A 0 A 1 A r 


Then like terms from equations (3.4) and (3.7) can be equated as follows. 

S a (8 ) = P 21 (A(8)P 11 ) i A(6)P 12 , i = 0, 1, ... , r 
a i 


(3.7) 


(3.8) 


The uncertainty modeling problem therefore requires that equations (3.8) be solved for P>i, P 12 , 

Pi 1 , and A( 8 ) such that the nilpotency condition of equation (3.3) is satisfied. 

In order to evaluate equations (3.8) and (3.3) in more detail, consider an expanded 

definition of P n , P 12 , and P 2 i containing partitioned submatrices associated with the 8 jl ni blocks 
of the A matrix given in equation (2.4a), as shown below. 
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P 12 = 


r 12 s 


(3.10) 


r 12 s 
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(3.11) 



where: 


r ll 


6 ; 6 ; 


GR 1 J . 


r 12s 


GR 


’cols 


r 21s 


gi j^ n rows xn i 


Equation (2.4a) is repeated here for convenience. 


(3.12) 


A = 


¥ni 



8 m I n r 


(3.13) 


Substituting equations (3.9) - (3.13) into equations (3.8) and (3.3) leads to a set of extremely 
complicated equations to solve. In order to satisfy the nilpotency condition of equation (3.3), the 
matrix P n must itself be nilpotent. Allowing P n to have a pre-defined nilpotent structure provides 
a means of somewhat simplifying these equations while assisting in satisfying the nilpotency 
condition of equation (3.3). The following Lemma establishes a general nilpotency structure that 
will be used throughout this paper. 

Lemma 3. 1 

Let A GR nxn be a quasi -triangular partitioned matrix whose main-diagonal blocks are nilpotent, 
as defined below. 


A = 


A 11 A 12 


0 A 


22 


0 0 

n ; 


Mm 

l 2m 


mm 


A M GR 1 *, Ajj = 0, rij <s nj , i = 1 , 2, ... , m 


Then matrix A is a nilpotent matrix with index of nilpotency, r|, as defined below. 

n m 

A = 0, r| = 2 "Hi ^ n 
i=l 


(3.14a) 


(3.14b) 


(3.15) 


Proof : 

Nilpotency of matrix A is clearly established by considering the eigenvalues of A. Since A is 
upper triangular, its eigenvalues are comprised of the eigenvalues of its main-diagonal blocks. 
Since each main-diagonal block is itself nilpotent, the eigenvalues of each must be zero (see 
Reference [7]). Hence, the eigenavalues of A must be zero and A must therefore be nilpotent. The 

index of nilpotency, r|, of matrix A is established by the following. 


Let: r = rj !+%+. + TJr 


A = A 


T ll + T l 2 + •••+ T l n 


n, n 9 r 

= A A •— A 
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nj 

Then, each matrix A contains a zero diagonal block corresponding to A n , since rj, is its index of 

j* 

nilpotency. It can therefore be shown that multiplication of these matrices to obatin A for r = rjj + 
r| 2 + ... + r| m results in the zero matrix, since each main-diagonal block is zero. However, if r < 
rii + ti 2 + ... + T) m then one of the main-diagonal blocks will not be zero, hence A will not equal 

zero. Thus, the nilpotency index for A must be equal to r = r| j + r\ 2 + ... + r) m . As can be 
verified in Reference [8], the nilpotency index for any matrix must be less than or equal to its 
dimension (i.e., n for matrix A). This can also be verified by the following. 

r|| <s n; for every i = 1, 2, ... , m 
m m 

=> *1 = ^Tii s ^lij = n 

i=l i=l 

Thus, equation (3.15) is satisfied. 

QED 

Note that the quasi-triangular structure defined by Lemma 3.1 is sufficient but not necessary for 
nilpotency. Other special structures can also be found. In fact, nilpotent matrices can be fully 

populated with nonzero elements. However, assuming some special structure for Pj j simplifies 
the solution of equations (3.8) and (3.3). For implementation purposes, allowing the special 
structure to be more general than upper-quasi-triangular may result in a less conservative (i.e., 
lower order) P-A model for some problems. However, for purposes of this paper, Lemma 3.1 

will be used to fix the structure of P, , so that the solution can be clearly derived. 

The quasi-triangular structure defined by Lemma 3.1 can be used in expanding equations 
(3.8) and (3.3). Thus, let P, , be defined to have the following upper quasi-triangular structure. 


where: 


p n, 

8 1 8 1 

Pi v 2 

P 1 L ft 
6 l 6 m 

0 

Pll8 2 8 2 ’ 

' ' P 1 L ft 

6 2 5 m 


0 

0 

‘ ‘ P 1 L ft 

6 m 6 m 


T >i 

(P. ) =0, Tli^nj , i=l,2, ... ,m 

ll8 i 8 i 


(3.16) 


(3.17) 


Then substitution of equations (3.10), (3.11) and (3.16) into equations (3.8) yields the following 
set of equations. 


Linear Terms: 



i = 1, 2, ... , m 


(3.18) 
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i = 1, 2, ... , m 


(3.19) 



Degree Terms: 







Crossterms: 


li -l li -l 

P 21 ( P ll ) 1 P 11 ( P 11 ) 2 ••• P 11 

21 S; 1 o; 8; n 6: 5: n 8; 5: 11 

*1 *1 *1 *1 *2 l 2 *2 


( p 1 


8 : 8j n 6i 6; 

nj -1 n r p nj n r p 


nj 


r 12 


8; 


ii'j’ 


= s 


) *(& ) 2 -(8: ) 

*1 *2 n T 


li 


(3.20) 


where: 


s =ii +ii +...+ii 


1 2 


h = 1, 2, ... , m - (n x - 1) 
i 2 = h + 1, ii + 2, ... , m - (n x -2) 


\ = ii + ( n T - 0, - > m 


n T = number of parameters in the crossterm 


Note that the S A terms on the right-hand side of equations (3.1 8) - (3.20) are the known constant 
matrix coefficients associated with the indicated parameter terms in S A (5). Moreover, depending 
on the number of parameters and the degree of each appearing in S A (5), there can be literally 
hundreds of S A coefficient terms - and hence equations to be solved. 


3.2 Numerical LFT Model Solution 

This section presents a numerical approach for solving all equations of the form defined by 
equations (3. 18) - (3.20) such that the nilpotency condition of equation (3.3) is satisfied and the 
resulting P-A model is of low-order. The results of this section are divided into three sub-sections. 
The first sub-section presents a solution for P 21 , P 12 , and the main-diagonal blocks of P u ; the 
second sub-section presents a solution for the off-diagonal blocks of P n ; and the third sub-section 
presents results relating to nilpotency and reducibility of the resulting model. 

3.2.1 Simultaneous Solution of P 21 , P I2 , and P u Main-Diagonal Blocks for each 5j Parameter 

The P 21 , P 12 , and P u main-diagonal blocks are solved simultaneously for each uncertain 

parameter 6; using the linear and § -degree terms defined by equations (3.18) and (3.19). 
Moreover, the solution is accomplished such that the resulting main-diagonal blocks of P, , are 
nilpotent with the appropriate index of nilpotency - as required by equation (3.17). This solution is 
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accomplished numerically with a matrix singular value decomposition (svd) by recognizing that 
this part of the problem is equivalent to a 1-D state-space (minimal) realization problem and by 

appropriately defining an equivalent Hankel matrix. The solution is accomplished for each 6; 
parameter as shown by the following theorem (which is based on Theorem 6-4, pages 268 - 272, 
of reference [9]). 

Theorem 3.1 


Consider the linear and ^-degree terms of Sa( 6) £ <p n rowsxn co b caQ ^ eX p anc j ec j as 

follows 


S V <8) = [ V + [S *i Si 2 ] 8 * + - + I s W 1 6 '’ 1 ' 


m 

,,= S [S A 

L >? n =l fl n-ls.n 


]5j n 


(3.21a) 

(3.21b) 


and use the constant coefficient matrices of equation (3.21) to form the Hankel matrices defined 
below 
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111 

u i 
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hj-1 
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hi 
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Using a matrix svd, factor equation (3.22) as follows 




o 

o 


0 0 
0 0 



= u 


ft. 2ft. V 6 . T = (U 8 . 2*. 1/2 )(Z X . 1/2 V x . t ) = P- 


a 'Si 


2U 12 


8; 


(3.22) 


(3.23) 


(3.24) 
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where: 


rank(S A ) = rankfP,, ) = rank(P 19 ) 

°6- 8 i 8 i 

Then the matrices P 21 , P 12 ,andP u form an irreducible realization of S A (6) as definedby 

5; 8; 8; 8; L/C 

equation (3.21), where: 


p^, =[i n oK, 

21 8j [ 'rows J 21 5j 

(3.25) 

ri n 1 

p = p n cols 

12 &i 12 8j 0 

(3.26) 

P 1L, = ( P 21, ^ S Aj ( P 12, ^ 

8 i 8 i 8; ig. 8j 

(3.27) 

and the notation (A/ designates the pseudoinverse of matrix A. 


Proof: 


From equation (3.24), define the following: 


<S. )' = (P, 2 ) f (P 2 l >’ 

% *>i *i 

(3.28) 

Then it is easy to show that: 


S* ^An ^An = P 21 P 12 ( P 12 Ml ) +P 21 P 12 

a 0 6 Z1 8j 1 Oj lz 8j 5j 8j 5j 

i 

= P71 Pi 2 = S* 

2 0; 12 8j A 0 6 . 

(3.29) 

Define the following relationship between the Hankel matrices of equations (3.22) and (3.23): 

S Al =m 8 .s a =s n 

§i 1 8j 6j 1 

which generalizes to: 

(3.30) 

V^o V ; n = 0 ’ 1 ’ 2 ’" 

8; 8; 

(3.31) 


where: 
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*rows 
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N 6; = 
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I n 0 0 0 
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0 I 0 

"cols 

; ; o i 

0 0 I 0 

"cols 


Consider the following: 

( P 1K ) 2 -[<P21 )%, (P]2 > t l 2 -< P 21 >%, < P 12 , Ml >%, <P,2. 

8 i 8 i 8 i l b: 8 i 8 i 1 6; 8 i 8 i l b: 


_ ( P 2U ^ M 6i S A 0i ( P 12 ft _ ^ M 6j S A 05 ( P 12 6 ^ 


' 5j ‘ 8 i 8 i 


8 i 


- ( P ?l )^M. Pi, P 19 (Pio )^(P?1 S* (Pit 


, '6; *21, *12, '■*12* > v 21, > ‘”6; A n '■*12* 

8 i 1 8 i 8 i 8 i 8 i 1 U 6j 8 i 


= (P 21 ) f M 6 . 2 S An (P ]2 ) f 


8 i 


8 i 


(P n ) n =(P 21 ) + M 6 . n S A (P 12 ) f 

6:6; 6; 8 1 A 0, JZ 6; 


Now, the constant coefficient matrices of equation (3.26) can be rewritten as follows: 


hi-1 


6; n 


I„ 0 

“rows 


IM 6i"-'SA 0 


6: 


11 cols 

0 


Substituting equation (3.29) into this expression yields: 


bi-1 


6; n 


I„ 0 

“rows 


M s . n - 1 S A „ (S. )'S 


*6; A 


0, 


n cols 

0 


Substitution of equation (3.31) into this equation yields the following: 


(3.32) 


(3.33) 


(3.34) 


(3.35) 
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n 


*n-l 


8; 11 


‘rows 


°k 0 vm 0 ) f s Ao ■"«* 

J 5: 8: 6: U 


Substituting equation (3.29) into this expression yields: 


Mi-1 


8; 11 




n rows J A 0 C v A 0 


An ( S An ^A* V^A* )% 


*0, ' A 0, 


11 cols 

0 


Substitution of equation (3.31) into this equation yields the following: 


k n-l 


8; 11 


rows 


•K 0 < S A 0 ^ )'S 

J u x. 1 u s u 5 "8; L 


8j 8 


I 


n cols 

0 


Substituting equations (3.24) and (3.28) into this equation yields the following result: 


S 


*n-l 


8j n 



rows 






Ml, > ,P 21 S 

8; 0: 


r 12 , 


k n-l 


8; n 


0 


‘rows 


P 21, « P 2I, ) tM 6i n_ 'S 4o < P I2„ )’> P ,2 

6 i 8j ' 


i “0, -8; ' ^8; 


n cols 

0 


Then, using equations (3.25) - (3.27) and (3.34) yields the following result: 


^ - P 21 (P[| ) P|2 > n — 2, ... , T]j 

n -l n 8j 8j8j 8j 


(3.36) 


Recalling equations (3.18) and (3.19), equation (3.36) shows that equations (3.25) - (3.27) are a 
realization of S A (6), as definedby equation (3.21). To show irreducibility, consider the 

L--r 

following: 


n j = dim(P 11§ § ) = rank(S A() 


min{rank(P 21 ),rank(P 12 )} 
8 : 8 : 


(3.37) 


Using equations (3.18) and (3.19), the following matrices can be defined to be consistent with the 
Hankel matrix given by equation (3.22) and its svd given by (3.24). 


n cols 

0 


11 



r 21 , 


21 Si 

(P n 8i8i ’ 

( V 8i > 2 


P2 > 8 < P " 8A ) ” i "‘ 


£i j^^1i n rows^ xn i 

6; 


(3.38) 


12 < 


f 12. (I'll.. > P 12, (I'll. . > 2p ,2. - ( p il. . )” i " lp i2. 

O; 0:0; O: 0:0: O: 0:0: O: 


p i2, eR 

o: 


n i x ^i n cols) 


(3.39) 


Since P 9 | GR*' 1 ' n rows ,xn ' and P 12 £R n ‘ x ^ llllcols) are tall and wide matrices (respectively) 

6j 8j 

that result from the svd computation of equation (3.24), the rank of each equals q and equation 
(3.37) can be evaluated as follows. 


n i =dim(P n ) = rank(S A ) = rank(P 21 )=rank(P ]2 ) 


V. 


Hence, the realization given by equations (3.25) - (3.27) is irreducible. 


(3.40) 

QED 


Note that as stated in equation (3.17), each main diagonal block of P, , must be nilpotent of 
index r];, i.e.: 

(P n ^ =0 

The following theorem establishes the nilpotency of P n 

8j8i 

Theorem 3.2 

The P n matrix computed using the result of Theorem 3.1 is nilpotent with index T|j. 

8 ; 8 ; 

Proof : 

Consider the following equation: 

<p ^ 8i ^ v t],,i 


12 



Substituting from equations (3.24) and (3.30), and using the fact thatUg and V 8 are unitary 

i i 

matrices yields the following result. 


('luj' 1 -[< U S:V' Z > ,M 8; S A„ (W'l 


1/2 yf 


1/2 V Vhhi 




<P > l» i6i , ’ li - ''X, X, V»i V »i V V“ ] ' 

fP dh _ ry -l/ 2 |j Tiy| it y y - 1 / 2 n^i 


Ty - 1/2 flj 


8 i 8 i 
1 6;S; 


( p iu , ^ =[2 6; - 1/2 U 6; Tm U.s, 172 ]' 11 


Then, the right-hand side of the equation can be separated into the product of matrix components as 
follows. 

(P (hi _ ry - 1 / 2 tj Tii| Ij y l/ 2 l 2 ry -l/ 2 |j Tiyj tt y 1/2 ihj “ 2 

^ rn 6-6- _L U Si iV1 5j 6j 6j J L ^6j u 5j iV1 5j 8 i 8 i 1 

Squaring the first term yields the following. 

(p llsi8i r =[( 2 8 i -nyM 5 u 8 z 8 / /2 )^ 

=* ( P IU . ) T1i =P 6i - 1 X i X i 2U 6 i S 6 i 1/2 HS 6 / 1 X i X i ^6 i S 6 i 1/2 ] T1i_2 


8 i 8 i 


Continuing this process yields the following result (which is consistent with equation (3.34) for n 
= *li)- 

(P,, ) T,i = 2. - 1/2 U fi T M ft hiU ft 1/2 

v 6; 5j 5; 6j 8j 


Since M a has r) i block rows and columns and is defined by equation (3.32), it is a nilpotent matrix 
i 

with index r\ { (see Reference [10]). Therefore, the deired result is obtained, i.e.: 

'Vi 

QED 


(P,,. . ) 1 =0 


In summary, this section has presented a simple numerical technique for computing P 21 , 
P 12 , and P u for each uncertain parameter. The result is irreducible, and each main-diagonal 

5; Sj 6; 

block is guaranteed to be nilpotent of index r|; , where r|j is the highest degree of 5j appearing in 
S a (6). 
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3.2.2 Solution of P n Off-Diagonal Blocks 


The P u off-diagonal blocks are each solved using the appropriate crossterms of S A (6), as 
defined by equation (3.20). The number of off-diagonal blocks to be solved is given by the 
following equation. 

m-1 

n ODB = (3.41) 

i=l 


The equation to be solved for each off-diagonal block of P, , is a generalized linear matrix 
equation. The general equation is given below for computing the off-diagonal block P u , 

8nSj 

where n = 1, 2, ... , m-1 and j = n+1, n+2, ... , m. 


< P 21. WP 1U M > P 11, . < P 12, )“S 4 

o„ o n dj §j o r 


[n] 


(3.42) 


n 


The matrices P, , M ? p M, p j2 , and S A M in equation (3.42) are comprised 
8 n 8 n 8 j 8 n 

of known matrices as well as matrices that have already been computed at this point in the solution 
process. Their explicit general definition is given in the following pages. 

The matrix P 21 ^ in equation (3.58) is a block-diagonal matrix with n partitions along 

8 n 

the main-diagonal, which is comprised of known matrices (i.e., matrices that have already been 
computed at this point). This matrix can be defined as follows. 


P 21 M = diag[P 21 ,P 21 [2] ,P 21 [2] ,P 21 [2] ,-, 

8 n 8 n 8 i j 8 n 8 ij 8 i 2 8 n 8 ij 8 i 2 8 i 3 8 n 


where: 


p [ 2 ] ... p 

r 21, o * , ’ r 21 

S il S i2 S ik-l S " 


8 1 8 2 " 8 i n . 1 6 n 


[ 2 ]] 


Partition 1: 


21 . 


21 s 

8 n 

P 21 s ( P lls s ) 

8 n 8 n 8 n 

P 21. ( p il. , > 2 

8 n 8 n 8 n 


( p in * ) 

8 n 8 n 8 n 


hn - 1 


n - 1 
0 


= 1 Block 


Partition 2: 


21 


6 i 1 6 n 


( 2 ( = diag[P 21 P],P 21 _ _ ra,.,P 2K _ PI] , 


21 s S ’* 21 , , 
8 l 8 n 8 2 8 n 


8 n-l 8 n 
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Partition 3: 


Partition 4: 



Partition k: 


n - 1 


= n-l Blocks 


P 91 [2] = P^j ® I ; ij = 1, 2, ... ,n-l 

Vn 8 ii 


r 21 


8 i 1 8 i 2 8 n 


|2| =diag|P 21 [2] ,P 21 [2] ,-,P 21 

8 l 8 2 8 n 21 6 1 5 3 5 n ^°l 8 n-l 8 n 


[2] 


[ 2 ] p [ 2 ] ... p 

1 9*71 9 9*21 

8 2 8 3 8 n 8 2 8 4 8 n 8 2 8 n-l 8 n 


[ 2 ] ... p 

’ ,r 2] 


8 n-2 6 n-l 8 n 


PI]; 


n-l\ (n - l)(n - 2) 

2! 


Blocks 


121 -p,, ®i„. „ 


21„ x x 21 x w ‘ 1 li^'4 n 

8 i 1 8 i 2 8 n 8 ij "2 111 


ij = 1, 2, ... , n-2 , i 2 = ii + 1, ij +2, , n - 1 


i 8 i 2 V" 


[ 2 ]=diag[P 21 [2] ,P 21 [2] ,-,P 21 

8 l 8 2 8 3 8 n 8 l 8 2 8 4 8 n 8 1 8 2 8 „-1 8 „ 


[ 2 ] 


[ 2 ] p [ 2 ] ... p [ 2 ] ... p 

9*21 ’ 9*21 9 9*2 


L 21 ~ ’*21 ~ ’"'’*21 ’"’’*21 1 » 
8 l 8 3 8 4 8 n 8 l 8 3 8 5 8 n 8 l 8 3 8 n-l 8 n “ 8 n-3 8 n-2 8 n-l 8 n 


11 “ (n - l)(n - 2)(n - 3) 
3! 


Blocks 


P I 2 1 = P (X) I 

21 S- 6- 6- 6 21 6- U \ 2 ^n ’ 

*1 *2 *3 n '1 z 


ij = 1, 2, ... ,n-3 ; i 2 = i l + 1, ij + 2, ... , n - 2 ; i 3 = i t + 2, i t + 3, ... , n - 1 
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[ 2 ] = 


21 


8 il 6 i 2 • 6 ik-l 6n 


diag[P 21 [2] ,P 21 [2 ] ,-,P 2 i 

8 l 8 2'” 8 k-2 8 k-l 8 n 8 l 8 2'" 8 k-2 8 k 8 n 8 l 8 2'" 8 k-2 8 n-l 8 n 


[ 2 ] 


r 21 


rap 


21 


[2] ... p 


21 < 


[ 2 ] 


8 1 8 3 ' 8 k-l 8 k 8 n 8 l 8 3’" 8 k-l 8 k+l 8 n 8 l 8 3”’ 8 k-l 8 n-l 8 n 


,P 


21 


8 n-3 8 n-2 8 n-l 8 n 


[2] ] ; 


n - 1\ _ (n — 1)! _ (n-l)(n-2)---(n~k + l) 

k-lj " (k - l)!(n - k)! ~ (k — 1)! 


Blocks 


E*21 ^ - P 21 • «• ri » 

5 ^ '\\2 • 'In 


Vi2 " Vl 6 " 


•l 


i, = 1, 2, ... , n - k + 1 ; i 2 = i, + 1, i, + 2, ... , n - k + 2 ; 


i k -2 = ii +k-3, ... , n -2 ; i = q + k - 2, ... , n - 1 


Partition n: 


21 


[ 2 ] _ 


= P 


8 l 8 2‘ 8 n- 


■l 8 n 


21 . 


01 


^2^3- • T l n -l T ln 



= 1 Block 


Note that all P 2U terms in the above equations are defined by the P 2U equation given for 
8 i 8 n 

Partition 1. 

The matrix P, , M in equation (3.42) is a block-column matrix with n partitions, and is 

8 n 

comprised of known matrices (i.e., matrices that have already been computed at this point). This 
matrix is defined as follows. 
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Partition 3: 


r ll 


[3] _ 


6 i 1 6 i 2 6 n 


r ll 


L 11 


[3] 


8 l 8 2 8 n 

I 

8 l 8 3 8 n 


[3] 


r ll 


[3] 


8 l 8 n-l 8 n 


r ll 


l ll 


[3] 


8 2 8 3 8 n 

I 

8 2 8 4 8 n 


[3] 


r ll 


[3] 


8 2 8 n-l 8 n 


r ll 


[3] 


8 n-2 8 n-l 8 n 


11 “ I'i (n-lXn-2) 
2! 


Blocks 


P„ [3] =P„ [2] P n [11 P n 

8 i 1 8 i 2 8 n 8 i i 8 i 2 8 n 8 i 2 8 n 8 n 


Pl^R, [21® I 


1 1 8 i 1 8 i 2 8 n " ll8 i 1 8 i 2 


p [ 2 ] _ p [l]p 

r l 1 s . "ML s r l 1 s 

V*2 S i, S i 2 6 i 2 


P„ 111 -P„ ®I„. 

6 ij*i 2 MM ' 2 


P U ' ' ' = P, , ®I„ 

8 i 2 8 n 8 i 2 8 n 1n 


ij = 1, 2, ... , n-2 ; i 2 = ij + 1, + 2, ... , n - 1 
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Partition 4: 


8 l 8 2 8 3 8 n 

I 

^ 8 l 8 2 8 4 8 n 


8 il 8 i2 8 i 3 8n 


8 l 8 2 8 n-l 8 n 

P 

ll8 l 8 3 8 4 8 n 

V 

ll8 l 8 3 8 5 8 n 


’ 3 


(n - l)(n - 2)(n - 3) 


8 l 8 3 8 n-l 8 n 


8 n-3 8 n-2 8 n-l 8 n 


Pit [4] =Pn [3] P n [1]P 11 

8 il 8 i 2 8i 3 8n 8 i 1 8 i 2 8 i 3 8 n 8 i 3 8 n 8 n 


P„ [3] =P„ [3] ®h 


p [3] _ p [2]p [ 1 lp 

r 1 1 — *1 1 r ] i *i l 

6:6:6: 1A S: 8 i^ 8 :^ 5: 5: “6s, 

*1 *2 *3 '1 '2 '3 *2 '3 *3 


etc. (see above) 


h = 1, 2, ... , n -3 ; i 2 = i] + 1, ij + 2, ... , n -2 ; i 3 = i t + 2, ij + 3, ... , n - 1 
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Partition k: 


8 l 8 2”- 8 k-2 8 k-l 8 n 

[1 

l8 l 8 2-" 6 k-2 8 k 6 n 


8 l 8 2” 8 k-2 8 n-l 8 n 


MS S ifc-l 6n 


8 1 8 3‘ " 8 k-l 8 k 8 n 


8 l 8 3'" 8 k-l 8 k+l 8 n 


8 l 8 3' 8 k-l 8 n-l 8 n 


8 n-k+l 8 n-k+2‘ " 8 n-l 8 n 


(n - 1)! (n-l)(n-2)—(n-k + l) 


k-i; (k-l)l(n-k)! 


(k-1)! 


Blocks 


. [k] =Pll lk - ,] Pll [1] P„ 

8 il 8 i 2 - 6 ik-l 6n 6 il 6 i 2 • 6 ik-l 6 n 11 «i k . 1 »n 8)1 


p n [kA] = Pu 

8 il 8 i 2 ^k-l 


p [k-1] _ p [k-2]p [l]p 

r l 1 “ r 11 r ] 1 r ] 1 

lip 5 p lip p p lip p lip 

O; O; ••*0; O; O; •••O; O; O; O; 

'1 x 2 ‘k-1 '1 >2 1 k-1 >k-2 1 k-1 *k-l 


P [ k- 2] =p [ k ‘ 2 ]®I 

11 6: 6: ■■•6; H 5| 5: -5: ^k-l 

‘1 l 2 ‘k-1 >1 >2 >k-2 


6:6:6: 
*1 *2 ‘3 


ra-P„ [2] ®I, 


p [2] _ p [l]p 

r ll s ~ “Ml. . r l 1 . 

6 : 0 : 0 : 0 : 6 : 

M '2 M *2 *2 
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p n [1] = p n ®K 

\\ '\% ”'2 


ij = 1, 2, ... , n — k + 1 ; i 2 = ij + 1, ij + 2, ... , n - k +2 


ik-2 = ii + k - 3, ... , n - 2 ; i k _j = i l + k - 2, ... , n - 1 


Partition n: 


P 1 1 M = p Pip 

8 l 8 2' 8 n-l 8 n 8 l 8 2-' 8 n-l 8 n 8 n-l 8 n 8 n 


^ ^ = 1 Block 


The first two matrices on the right side of the above equation for Partition n are defined by the 
preceeding equations for Partition k. Also, all P, , terms in the above equations are defined by 

8 i 

the P| | equation given for Partition 2. 

8 n 

The matrix P 12 in equation (3.42) is a block-row matrix with j partitions, and is 
8 j 

comprised of known matrices (i.e., matrices that have already been computed at this point). This 
matrix can be defined as follows. 


12 . 


“L P 12- P 1 


11 


Sj 8j 


12 , 


’( P l 


1 ) 2p 

¥) 


12 , 


>( P 1 


1 ) 11j " 

¥) 


r 12 , 


The matrix S A ^ 


on the right side of equation (3.42) is a block-column matrix with n 


partitions, and is comprised of known coefficeint matrices from the expansion of S A (6). This 
matrix can be defined as follows. 
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s A [11] = 


X 8 j 

s A pi 
Nvj 
s A ra 

XjVn 8 ] 

s A P] 


*1 >2 >3 n J 


[ 2 ] 


V'2"V n 4 j 


[ 2 ] 


8 l 8 2"' 8 n-l 8 n 8 j 


Partition 1: 


Vj 


8 n 8 j 


V 8 j 


2 x a 2 
6 n 8 j 


V 8 j 2 


SA % SA (^n +1 ) 

8 n T ln8j 8 n r in8j^ 


8 n 8 j 11 J 


^(T|j +1) 


6 n 26 j T| j 


Hrin+rij-!) 


6 n T|n8 j r| j 


Partition 2: 


j = n + 1, n + 2, , m 


[2] _ 




\Wj 

\*n 2 *j 


8 i 1 8 n Tln8 j 


; ij = 1, 2, ... , n - 1 ; j = n + 1, n + 2, 
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VA 


\i+0 


6 ;i 6./6j 


H2+t) 2 

Vn'aj 


v2+t?) 




v (3+/) 


6- 2 6 ^S- 2 
u ij u n u j 


A (ilij +0 


Sij^lSn^j 


A (rij i +e+l) 


5 i 1 11i l 5 n^ 6 * j 2 


Partition 3: 


THj+f) 


v„v 


'J 


>A (Tlj+€+l) 


JijWi 


>A (rii +rj j +i-\) 


8 i 1 11i l 8 nV 1j 


; ^ = 1,2, ... , T| r 


[ 2 ] _ 


8 ij 8 i 2 8 n 8 j 


S i, S i 2 < ’ i 2Vj 

is i 1 S i 2 fi 2 6 „ 2& j 


6 i 1 6 i 2 < ’ i 2 6 n’'" 6 j 


ij = 1,2, ... ,n-2 ; i 2 = i t + 1, ij + 2, ... , n - 1 

4 = 1, 2, ... , r| j ; j = n + 1, n + 2, ... , m 

2 2 
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’A - 

6^ 


J> 

+ 

^ 1 

1 

A (2+f) _ 2 

13 A 

Olj+f) 

6 

8 i i 88 j 

8 i i 88 j 


S A - 

- S A, - 

(2+*) 

V“j 


(ilj + £+l) 

$ 


>1 


(ilij +€) (r^+f+1) 

V^i 


*1,^1 “j 2 


(rijj +rij+£-l) 


Partition 4: 


8 = ( 8 iJ i2 (\/ n ; K = 1 * 2, ... , T| n ; £-£. 


Note : £\ is updated before £ n 


[ 2 ] _ 


8 i 1 8 i 2 8 i 3 8 n 8 j 


8 i 1 8 i 2 ^ i 2 8 i 3 £ ‘3 8 n 8 j 
^ 8 i 1 8 i 2 £i 2 8 i 3 ^ i 3 8 n 28 j 


8 i 1 8 i 2 £ '2 8 i 3 fi 3 8 n t1n8 j 
ij = 1,2, ... ,n-3 ; i 2 = i t + 1, ij + 2, ... , n - 2 
i 3 = ij + 2, i l + 3, ... , n - 1 ; j = n + 1, n + 2, , m 


4 = 1, 2, ... , ri j ; £ { = 1, 2, ... , r| j 

2 2 3 3 


8 VJ 

, 288 j^j 

i^l 88 ^ 
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7 A - 

V 8 J 


V<> - 

»i,“j 

- 

(2+0 - 2 
V 8 ] 

'(2+1) 

S A - 


Olj+O 


8 i 1 88 j 11 j 


(rij + ^+1) 


8 ij 288 j^j 


(Tlij +0 


V , .. -Sa 




(ri^+f+1) (riij +r|j + £-l) 


Partition k: 


6 = (6 i2 ) 12 (6 i3 ) ‘ 3 (6 n / n ; £ n = l, 2, ... , r\ n ; *-4 2 +4 3 +* n 


Note : £ \ is updated before ^ ^ is updated before 


[ 2 ] _ 


V‘2 "Vi 8 " 8 ; 


6 il 6 i 2 (i 2 ' ^k./ik-lVj 

'i 

A 6i i 5i 2 ^ 2 "- 6 i k _ ] \-i 8 n 28 j 


8 i 1 8 i 2 £i 2 - 8 i k ./'k-l 8 n T1n8 j 

ij = 1, 2, ... , n - £ + 1 ; i 2 = ij + 1, ij + 2, ... , n - £ + 2 


k -2 = ii + - 3, ij + ^ - 2, ... , n - 2 ; y j = ij + £ - 2, ij + ^ - 1, ... , n - 1 

j = n + 1, n + 2, ... , m 


£\ = 1, 2, ... , ri j ; ... ; =1, 2, , r| j 

2 2 k-i k-i 
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’A - 

6^ 


V<) - 
V 8 * * j 

A (2+f) _ 2 

V 6 J 

'(2.7) 

V“j 

S A - 

<3 *\^ 


Olj+f) 




(rij + £+l) 


6ij 2 66j T 1j 


Olij +£) 


V -S 4 


Vl“j 


(r^+f+l) (riij +T lj + ^-l) 

Vl“j 2 


8-(»i 2 ) ' 2 -(8 ik _ 1 ) ik -'(8„)' n ; <„ = 1. 2 «1„ 


£ = £-. + £■. + ••• + £■ + £ 

l 2 >3 'k— 1 11 


Note : L is updated before L ; L is updated before £. ; ... ; L is updated before £ n 

2 3 3 4 k-1 


Partition n: 


[ 2 ] _ 


6 1 8 2 " 6 n-l 6 n 6 j 


8 1 8 2^ 2 "' 8 n-r n-18 n 8 j 
i 5 1 5 2 ^2---5 n _/n-l5 n 2 5 j 


8 1 8 2^ 2 ‘" 8 n-l^ n_l8 n^ n8 j 

j = n + 1, n + 2, ... , m 

£ 2 = 1, 2, ... , T] 2 ; ••• l £ n -l = 1, 2, ... , T| n _! 
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’A - 

6j66j 


+ 

^1 

1 

- 

(2+0 - 2 

6j66j 

6j66j 


in 

> 

1 

(2+0 

6j 2 66j 

V^ 2 


S A - s A _ 

Oll+O (ni + Ml) 

6^2 


Olj+O 




blj+£+l) 




(ril+rij+f- 1 ) 




6 = (6 2 ) 2 ---(6 n _ 1 ) n -l(6 n ) n ; £ n = 1,2,... ,ri n 

J = l 2 + h + "‘ + ^n-l + ^n 

Note : l 2 is updated before is updated before i 4 ; ... ; l n _\ is updated before £ n 

The above general equations, which define the matrices given in equations (3.42) for 
generating the off-diagonal block equations, are complicated due to the large number of cross- 
product terms that can arise in solving the general problem and due to the notation required to 
generate the associated equations. As an illustration of generating these equations based on 
equations (3.42) and the above defining equations, the off-diagonal block equations for the case of 

three parameters (m = 3) with maximum degree of 2 for each 6j parameter (r| { = r\ 2 - % = 2) are 
shown below. 


Off-Diagonal Block Equations for m = 3 (r| j = t] 2 = b ? = 2> 

S, 


r 21 , 


>1 


P 2 1 P 11 

“Oj Oj6j 


Tl 


8 l 8 j 


Pj2 P 11 P 12 

8 j ¥i 8 j 




' 8 l 8 j 2 


38 r 8 r 


, j = l,2 


21c 

P 2U P 11 


S 2 S 2 6 2 


P 21 c P 11 


6j “6 j6 2 


P 21c P llc S Mlc c 
O] 0l°l O]0 2 


P 1 
P 1 

P,, P, 


P 21 c P llc c *11 c c 

O) 0102 o 2 o 2 


P 21 c P llc c Ml* c Mlc c 
Oj ojOj 0j02 O90 


2°2 


Tl 


8 2 8 j 


P 12 P 11 P 12 
S ¥i *j 


S 4 . 

‘ S 2 6 J 

S A 1 

A 2 o 

8 2 8 j 2 

S A 

A 2 o 

8 2 “ 8 j 

Sa 3 

W 

Sa 2 

2 8 i 6 2 6 j 

% 2 
6 i 6 2 6 j 2 

S A 
a 3 o 

8 i 2 8 2 8 j 

S A 

a 4 0 0 

8 i 2 6 2 8 j 2 

Sa 3 

6 i 8 2 2 8 j 

Sa 4 

V2V 

S A 

V^j 

Sa 5 

8 i 2 8 2 2 8 j 2 
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The general equation for computing each Pj j off-diagonal block, P n , given by 

5 n 6j 

equation (3.42), can be written as a generalized linear matrix equation of the following form: 

AXB = C (3.43) 

where A, B, and C are known constant matrices. The following Lemma is stated without proof as 
an extension of Lemma 2.2 given in Reference [11]. 

Lemma 3.2 


Consider the generalized linear matrix equation given by equation (3.43), where A G R nxm , B G 
R“ p , and C G R nxp are given matrices. Then the following statements are equivalent: 


mxr 

(1) there exists a solution X G R ; 

T 

(2) the columns of C G Im [A] and the rows of C G Im [B ]; 

(3) rank [ A C ] = rank [ A ] and rank [ B T C T ] T = rank [ B ]; 

(4) Ker (A*) C Ker (C*) and Ker (B) C Ker (C). 

Furthermore, the solution, if it exists, is unique if and only if A has full column rank and B has 
full row rank. 


Equation (3.43) and Lemma 3.2 can be used in computing a solution for each off-diagonal block of 
P n , based on equation (3.42). This solution has the following form. 


where: 


4 


M = B ® A 

4 • 


M\N 


N =C 


1 


Note : C is the column-form vector of matrix C obtained by 

vector 

i = 1, 2, ... , r 


stacking the columns of C into one column vector 


X = [ x/ X 2 j ... X r l ] ; X; 1 G R mx 


(3.44) 

(3.45) 


(3.46) 


Then the following theorem is stated. 
Theorem 3.3 


Given a general linear matrix equation of the form given by equation (3.42) for each off-diagonal 
block of Pj j, i.e.: 


(P 21 MP„ M)P„ (P 12 )-S, 

6 n o n 6 n 6 j 


[n] 


where: n = 1, 2, ... , m-1 and j = n+1, n+2, ... , m 

then a solution for P n of the form given by equations (3.43) - (3.46) and which satisfies rank 

6n6j 

test (3) of Lemma 3.2 always exists and is irreducible. 
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Proof : (Sketch) 

The rank test (3) of Lemma 3.2 can be used to determine whether a solution for P n exists, 

n 6nSj 

based on the P 21 , P 12 , andP u matrices computed as described in Section 3.2.1. If not, 

6 ; 6 ; 5 ; 6 ; 

these matrices can be augmented using the appropriate columns and/or rows of the matrix S A ^ 

given on the right side of equation (3.42). Thus, a solution can always be found. The resulting 
solution is irreducible, because satisfaction of rank condition (3) in obtaining a solution prevents 
unnecessary redundancy from being built into the solution process. 

QED 


To summarize, this section has presented a simple numerical technique for computing the 
off-diagonal blocks of P n , i.e. P n , for each block-row, n, and each block-column, j, (as 

S n 6j 

defined by equation (3.16)), where n = 1, 2, ... , m-1 and j = n+1, n+2, ... , m. The numerical 
computation involves the solution of a generalized linear matrix equation, and such a solution can 
always be found by augmenting the previously computed P 2] , P 12 , P n , and P n 

5j 6j 6j 5j 5; 6j 

matrices as required to obtain a solution for equation (3.42) based on equations (3.43) - (3.46). 
The result is irreducible, because a solution for each off-diagonal block is computed to just meet 
the rank conditions (3) given by Lemma 3.2. 


3.2.3 Full P-A Model Solution, Nilpotency and 1-D Irreducibility 

Once the P 21 s * P 12* ’ P 1U s , and P n partitions for each parameter have been 

djdj 6idj 

determined as described in Sections 3.2.1 and 3.2.2, the full solution is determined using 
equations (3.9) - (3.12). This is a simple matter of collecting the matrix partitions together into a 
single matrix for P 21» P 12’ and Pjj. The A matrix is also known and given by equation (3.13), 

where the number of repetitions for each parameter, nj, was determined in solving the P 21, > P 12, > 

6; 5; 

Pi i , and Pi i matrices. 

8i5i 1 oj6j 

The following theorem is given regarding the satisfaction of the nilpotency condition of 
equation (3.3) for the full P-A model solution. 

Theorem 3.4 


The P n matrix defined by equation (3.9) and computed using Theorems 3.1 and 3.3 as described 
in Sections 3.2.1 and 3.2.2 satisfies the nilpotency condition of equation (3.3), as defined below. 

[ A Pjj ] r+1 = 0 ; r+1 < t)j + r| 2 + ... + rj m 


Proof : (Sketch) 

For r+1 = r|j + r| 2 + ... + r| m , nilpotency is satisfied by Lemma 3.1. For this case, solution of the 
off-diagonal blocks does not enter into satisfying the nilpotency condition. That is, the nilpotency 
of the main-diagonal blocks is sufficient to satisfy the nilpotency of the full solution. 
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For r+1 < r| 1 + r\ 2 + ... + r| m , nilpotency is satisfied by the solution of the off-diagonal blocks. 
That is, this case arises when there are zero crossterm coefficient matrices that are factored into the 
solution of the off-diagonal blocks. Thus, inclusion of these zero matrix coefficients in the 
solution of the off-diagonal blocks automatically satisfies the nilpotency of the full solution. 

QED 


An objective of the P-A modeling process was to determine a model which is low-order. 
The following theorem is therefore given regarding the reducibility of the full P-A model solution. 

Theorem 3.5 

The P-A model matrices defined by equations (3.9) - (3.13) and solved using Theorems 3.1 and 
3.3 is 1-D Irreducible. 

Proof : (Sketch) 

The P 21 , P 12 . , and P n matrices determined using Theorem 3.1 represent an irreducible 

5i 6i 5i5i 

realization of the linear and n -degree terms of S A (6) associated with the 6, parameter. Solving 
equation (3.42) using Theorem 3.3 results in an irreducible solution of the off-diagonal blocks of 
P u based on the solution obtained previously for P 21 s .> P 12, , and P H . Thus, putting the full 

8i 6i 5j5i 

solution together results in a 1-D irreducible LFT model of the given system. 

QED 


4. Example: Multivariate Quadratic Problem (See Reference [4]) 

Consider the following compound inertia matrix problem presented in |4], and first posed 

in [13]. 


0 


J = 


2yz 

-2y 2 


-2yz 

2y 2 

4(y 2 - z 2 

0 

-2xy 

-4xy 

2xy 

0 

4xz 


-3xy xz 

3(x 2 - z 2 ) yz 

a 2 2 
-3yz y - x 


(4.1) 


The x, y, and z terms represent displacement parameters from some reference (zero) point for the 

system. Thus, the parameters x, y, and z are the uncertain parameters, 5, of the system. The 
results obtained using the above computational solution (in Matlab) are shown below. However, 
the details of obtaining this solution are omitted for brevity. 


P 21 “ [ P 21, 1*21, P 21, ] 

Ox Oy Oz 


r 21 


Sx 


■ 0 0 0 0 1 ' 
1.7321 0 0 0 0, 
0 0 10 0 



0 2.1147 0 0 0 0 O' 
000 0 000 
0 0 0 -1.4953 0 0 0 


(4.2a) 


(4.2b) 


30 



r 21 < 


20 0 000 

0 0 1.7321 0 0 0 
00 0 001 


(4.2c) 


*12 = 


r 12 . 


0 0 0 0 0 0 

0 0 0 0 1.7321 0 

0 0 0 0 0 0 

0 0 0 0 0 1 

0 0 0 0 0 0 


r 12, 


r 12 , 


r 12 , 


12 < 


r 12 , 


(4.3a) 


0 
0 

1.3375 0 
0 0 
0 0 
0 1 
1 0 


0 .94574 1.8915 0 0 

0 0 0 0 0 


0 0 0 0 0 0 

0 0 0 2 0 0 

0 0 0 0 0 0 

0 0 0 0 1.7321 0 

0 0 0 0 0 1 

0 0 0 0 0 0 


0 

0 

0 

0 

0 


0 

0 

0 

0 

0 


0 -.66874 
0 0 

1 0 

0 0 
0 1 


(4.3b) 


(4.3c) 
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8 X 8 X 


*11 = 


* 11 , , * 11 , , *11 

O x O x o v o 


0 

0 


r ll 


'x~y 


8y8y 


8 x 8 z 

0 


*11 *11 

8 z 8 y 8 z 8 z 


(4.4a) 
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(4.4b) 
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r- 1.2209 0 

0 

0 

0 

0 

O' 


0 

0 

0 

0 

0 

0- 

0 

0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

0 

, Pll 

5 V 5, 

0 

2 

0 

0 

0 

0 

0 

0 

0 
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0 

0 

0 

X z 

0 
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0 

0 

0 

0 

0 

0 

0 

0 

-3 

0 

0 


0 

0 

0 

0 

1 

0 


r ll 


6 z 6y 


0 0 0 0 0 -1 0 ■ 

0 0 0 0 0 0 0 

0 0 .28778 0 0 0 .7698 

0 0 0 0 0 0 0 

0 0 0 0 0 0 0 

0 0 0 0 -3 0 0 


(4.4c) 


(4.4d) 


A = diag [6 X I 5 , S y I 7 , S Z I 6 ] (n A =18) 


(4.5) 


Note that the solution of this problem was not restricted to a quasi-upper-triangular Pj l matrix. In 
particular, it was determined in solving this problem that the quasi-upper-triangular structure for 
P) j required an extra repetition in A to obtain a solution. 

A comparison of this solution with those obtained in [4] and 1 13] is shown in Table 1 . 


Method 

EH 

eh 


n A 

Comments 

Belcastro & Chang 

5 

7 

6 

18 

Direct Numerical Solution for Nonlinear 
Problem, No Decomposition, No Model 
Reduction 

Cockbum & Morton [4] 

9 

10 

9 

28 

Decomposition to Linear Components, 
Solution for Each Linear Component, 
Combination of Component Solutions 


7 

8 

5 

20 

Same As Above with Model Reduction 


7 

5 

7 

19 

Special Decomp, to Linear Components, 
Solution for Each Linear Component, 
Combination of Component Solutions 


6 

5 

6 

17 

Same As Above with Model Reduction 

Doyle, Elgersma, et. al. [13] 

9 

9 

9 

27 

Decomposition to Linear Components, 
Solution for Each Linear Component, 
Combination of Component Solutions 


I 

5 

1 

13 

Special Matrix Decomp, to Linear 
Products, Solution for Each Linear 
Component, Combination of 
Component Solutions 


Table 1. Comparison of LFT Models Obtained Using Current Methods 
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The solution obtained using this LFT modeling approach required a total of 18 parameters in A, 
with 5 repetitions for 6 X (n x = 5), 7 for 6 y (n y = 7), and 6 for 6 Z (n z = 6). Note that the LFT 
modeling approach of this paper does not require matrix decompositions for a solution to this 
example, since it was already in a multivariate polynomial form. Moreover, this approach achieves 
a low-order model directly (without the use of model reduction), and can be readily implemented in 
Matlab. The result presented in [4] for a direct decomposition required 28 and 20 parameters in A 
before and after model resuction, respectively. The result obtained using a specialized 
decomposition approach developed in [4] to reduce the resulting LFT model dimension required 19 
and 17 parameters in A before and after model reduction, respectively. Note that this approach 
decomposed the J matrix of equation (4. 1) to linear matrix products and sums. Then an LFT 
model for each linear matrix was obtained separately, the individual LFT models combined to form 
the full LFT model, and reduction methods applied to remove unnecessary repetitions. The result 
presented in [13] required 27 parameters in A using a linear decomposition approach, and 13 
parameters in A by recognizing that J can be factoed into the product of two matrices containing 
only linear x, y, and z terms. Although this yields the lowest-order LFT model, it is specific for 
this particular matrix structure and can therefore not be generally applied to other problems. 


5. Concluding Remarks 

A numerical approach was presented in this paper to directly compute low-order LFT 
models for multivariate polynomial problems. The LFT modeling approach does not require 
matrix decompositions for multivariate polynomial problems, and a low-order model is directly 
obtained without model reduction. The computations depend only on simple matrix computations, 
including the singular value decomposition (svd) and solving generalized linear matrix equations. 

A matrix svd is used to simultaneously compute a solution for the P 21 , P 12 , and P n 

6; 6j Sj 6j 

matrices for each 6[ parameter. Generalized linear matrix equations are used to solve for the 

P n matrices. The full LFT model is constructed by simply collecting the partitioned solutions 

8;5j 

together into the P21’ Pl2> and P n matrices. The resulting LFT model is low-order, because 
matrix structure is exploited during the computations in satisfying the rank conditions required for 
a solution. Future work will include developing a Matlab implementation of this LFT modeling 
approach. 
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